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S ■ Abstract 

Baryogenesis at the electroweak phase transition may take place through CP- violating reflec- 
tions of quarks from expanding bubbles of the broken symmetry phase. We formulate and 
approximately solve the transport equations for the reflected quark asymmetries. The results 
are comparable to a previous Monte Carlo calculation and allow for analytic estimates of the 
baryon asymmetry AB when the bubble walls have a small velocity. Our method predicts 
no velocity-dependence for AB in the small v limit, unlike results obtained from a simplified 
treatment based on the diffusion equation. 



There is currently much interest in the possibility that the baryon asymmetry of the 
universe was created during the electroweak phase transition, using its first-order nature as 
the means for departure from thermal equilibrium, as is required in any theory of baryo- 
genesis [1-4]. There are two paradigms, depending on whether the bubbles of true vacuum 
which nucleate in the symmetric phase during the transition have relatively thin or thick 
walls. In the thin-wall regime which is investigated here, CP-violation in the Higgs sector 
causes an asymmetry between fermions and antifermions reflected from the bubble walls. 
To leading order in a wca k it consists of equal and opposite excesses of right-handed and left- 
handed fermions, due to CPT conservation, so that the net baryon or lepton number in the 
disturbance is zero, but its chirality is nonzero. This asymmetry is subsequently converted 
to a baryon asymmetry by sphaleron interactions. The expanding wall overtakes the baryon 
asymmetry and incorporates it into the true vacuum phase where it is preserved because of 
the freezeout of sphalerons. 

To compute the resulting baryon asymmetry in this scenario, it is necessary to know the 
distribution functions of the fermions and antifermions in front of the wall in some detail. So 
far several treatments have been given, including a Monte Carlo simulation of the diffusion 
of reflected quarks || (hereafter called CKN), and the diffusion equation M. Our aim is to 
complement these approaches by starting with the exact Boltzmann equation and solving for 
the distributions with the help of a few reasonable approximations. We will find the profile 
for the reflected quark asymmetry, estimate the effect of nonvanishing domain wall thickness 
on the baryon asymmetry, and show that the present method disagrees with the diffusion 
equation approach concerning the velocity dependence of the final result. 

Our starting point to compute the transport of reflected quarks with distribution func- 
tions fi is the Boltzmann equation, 

/i(x,p) = C, (1) 

where C, is the collision term.[] The latter is the standard difference between two terms, 

1 We have ignored force terms that would give rise to the screening of hypercharge || because any de- 
viations Sfi due to screening would be of the form yi<f>, where yt is the hypercharge of particle i and cp is 
the hypercharge potential to be screened. Such a deviation would have no effect on the biasing of sphaleron 
interactions because they conserve hypercharge. 
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which scatter the particle of interest respectively into and out of an element of phase space, 
Q = Yl [dn l dn 2 dnJ(27r)*8\p 1 +p 2 -ps- P )\Mi2^3 P \ 2 f{p2)fi(x,Pi) 

channels 1 

-(Pi<-p)}; d^ = d 3 Pj /(167T 2 E,), (2) 

ignoring Pauli blocking factors, /(p) is the nearly-equilibrium distribution function for the 
quark or gluon off of which the species i of interest is scattering. Of course there is one of 
approximately 20 channels where the former is also the quark of flavor i, in which case fi 
would appear quadratically, but ignoring this makes only a small error. We have also omitted 
decay processes because the thermal masses of the gauge bosons are below the threshold for 
producing two thermal quarks, and the rate for Higgs boson decays into fermions is smaller 
than scatterings mediated by gauge bosons, due to the large number of scattering channels. 

Because the scattering processes are dominated by low momentum transfer, due to the 
near masslessness of exchanged gauge bosons in the thermal plasma, it is possible to simplify 
the collision term using the Fokker-Planck approximation |J. The key observation is that 
in the absence of a thermal mass m g for the gluon, Cj would diverge logarithmically in the 
infrared because the matrix element has the leading behavior \Ai\ 2 ~ g^s 2 / '(t — m 2 ) in terms 
of Mandelstam variables. We will isolate the logarithm and discard the corrections that 
would vanish as m g — > 0. Let p denote the three-momentum of the scattered quark, and 
define P2 = p'; P3 = p' + q; so p' is the momentum of the background quark or gluon and 
q is the momentum transfer. By substituting p x = p + q in the first term on the right hand 
side of @j and pi = p — q in the second term, and doing the integral over c?n 3 using the 
momentum-conserving delta function, we get 

Ci = I ^p/rf 3 g{G(p + q,q)-G(p,q)}; (3) 

G( P ,q) = i ^ P) i / 1 d nf f }^ i\ gflpl + 1p 7 I - Ip - ql - 1p 7 + a\)\M p , p ^ p - q , p , +q \ 2 
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as the contribution from a single scattering channel; we will sum over channels in the end. 
The next step is to approximate the integrand of ([3D as 
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by Taylor-expanding. This is a good approximation since we expect the scattering to be 
dominated by small momentum transfer. One finds that only these first two terms are 
logarithmically sensitive to the gluon mass; the higher ones are convergent. It is farily 
straightforward to evaluate them, with the result that the collision term can finally be 
expressed as 

Q = Dd p ■ ^,/,(x, p) + £H/ 4 (x, p)) • (5) 

This form is quite general and appears in nonrelativistic scattering processes as well . 
The diffusion coefficient emerging from this reduction is 

D = »^ I (6) 

where the factor of 20 comes from weighting the multiplicity of different scattering channels 
by their respective squared matrix elements and thermal statistical factors (8x1x1 for 
gluons, 18 x 4/9 x 3/4 for quarks and the same for antiquarks). To evaluate the logarithm 
we take the Debye screening mass m 2 = 8ira s T 2 for the lower limit, the thermally averaged, 
squared, center of mass energy s = 20T 2 for the upper limit, and a s = 0.1, with the 
result D = 0.16T 3 . Since we have ignored weak interactions, D is the same for both quark 
chiralities. Note that D differs in dimension and meaning from the usual diffusion coefficient, 
since it multiplies derivatives with respect to momentum rather than position. 

It is convenient to work in the rest frame of the wall, where the distribution will be 
stationary so that the time derivative can be neglected. In this frame we should Lorentz- 
transform the momentum variables in (|l|) and (JD, but this complicates the solution of the 
transport equation, so we will take the limit of small wall velocities v, where the difference 
between the momenta in the rest frame of the universe and that of the bubble wall can be 
neglected. It will turn out that the asymmetry we seek between quarks and antiquarks is of 
order v so that simply dropping the time derivative in eq. ([!]) amounts to ignoring higher 
order corrections in v. 

The term linear in derivatives in eq. (|5|) is cancelled if we write fi = e~ l3E / 2 fi. Also, since 
we are interested in the asymmetries between particles and antiparticles, 5f = fa — fj, we 
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must subtract from eq. (ffl) the same equation for /j. The solution will have the form 

Sf(p, z) = e-^ 2 £ r dke~ k ^ (d(z)At n gt n (p) + B{-z)A^gJJp)) . (7) 

n=0 Jk ° 

The differential equation for gf n is 

(^~4^ + 2 \e~E^) T De) 9 ^ = ^ (8) 

where one must remember that E = p for g£ n but = (p 2 + m 2 ) 1 / 2 for because the 
electroweak symmetry is broken only to the left of the wall. However, we will solve (§) only 
to leading order in the quark mass m. Then, similarly to our neglect of the wall velocity in 
the Fokker-Planck equation, we can set E = p to discard higher order corrections in m, since 
the solution 8f is homogeneous in a source term that vanishes as m — > 0. Throughout this 
paper m is understood to be the quark mass in the broken symmetry phase at the critical 
temperature T c , which is smaller than the usual mass because the Higgs field VEV is smaller 
at T c than at T = 0. 

Using parabolic coordinates p z = |(£ — 77), p — +T)), the massless limit of eq. (||) is 
separable and it can be solved exactly: 

9kn(^v) = e- a ^ 2+ ^ 2 1 F 1 (-n;l;a + r ] ) 1 F 1 {c kn -l--ia^); (9) 
a± = (k/D±f3 2 /A)^ 2 ; c kn = 1/2 - (i/a_) (J3/2 - (n + $a+) , 

and gknUiV) = QkniViO- The requirement that n be an integer in eq. @ is to insure good 
behavior as r] — > 00. The functions g^ n form a complete set for k > k = (3 2 D/A = 0.04T. 
Note that the spatial extent of the diffusion layer will therefore be given by the distance 
scale fcg , from eq. (0). 

Let us now discuss the boundary conditions for the distributions, which express the 
conservation of flux at the bubble wall. For example, since left-handed particles are reflected 
into right-handed particles and vice versa (because of angular momentum conservation), the 
distribution function for left-handed quarks satisfies 

f+( Pz ) = n R ^ L f+(-p z ) + r L f-( Pz ) ( Pz > m), (10) 
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where the superscripts indicate that fi is evaluated at z — 0+ or z = 0— with respect to 
the wall, TZr~>l is the probability for the reflection of right- to left-handed quarks, and Tl 
is the transmission probability for left-handed quarks. Subtracting the distribution function 
of the antiparticle to obtain Sf + = f^ L — remembering that has the corresponding 
asymmetry —5f + , and doing the same for the distributions on the z < side of the wall, one 
obtains 

5f + (±p z ) - 5f-(± Pz ) = tK (5f + (- Pz ) + 5f-(p z )) + S(p z ) (p z > m), (11) 

where IZo is the mean value of the two reflection coefficients IZr^l, T^-l->r, and the source 
term S(p z ) depends on their difference A7Z, and on that of the unperturbed distribution 
functions: 



s(p M ) = (f-(p z )-f + (-p z ))An( Pz ) 

S 2{3v Pz e-P p An(p z ). (12) 

Here we used the equilibrium distribution functions f(±p z ) = e~ l3(j ' ±vpz ' > on either side of the 
wall, expanding to first order in the velocity v. The difference between positive and negative 
p z is due to our working in the rest frame of the bubble wall. 

As one would intuitively expect, the solution to the boundary conditions (|IT]) requires 
that Sf~(p z ) = —5f + (—p z ): there is a dipole layer of left-handed quark excess at the bubble 
wall, since the net number of chiral quarks is conserved and remains zero during the process 
of reflection. This means that A^ n = — A^ n , and the boundary condition becomes an integral 
equation for A^ n , which can be solved using the orthogonality of the eigenfunctions rj): 

Kn = N£ [°° f* drj^ - v)e M2 S(p(i, r,)) [g^t, rj) - g+ n ( V , £)] ■ (13) 
jo Jo 1 J 

Here N kn is the normalization of the orthogonality integral of the g^'s, N kn = (D/Sa + ) sin(7TCfc r , 

We thus have an approximate solution for the left-handed quark asymmetry in the sym- 
metric phase, which enables us to compute the baryon asymmetry that sphalerons will be 
driven to produce. A convenient measure is the baryon-to-entropy ratio. In CKN it is found 
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that in a theory with n Higgs doublets (we will assume n = 2 

10- 7 /€ 



(2n + l)u 



/^tfn, (14) 



where k is the dimensionless constant (k < 1) characterizing the rate of sphaleron interactions 
in the symmetric phase, and 6n is the number of reflected quarks per unit area in front of 
the wall, 



5n 



/roo 
d 3 p / dzSf(p,z). (15) 
J ZQ 



Notice that in our approximation of small wall velocities, 8n ~ S(p) ~ v so the baryon 
asymmetry is independent of v in this limit. 

The lower limit zq of the spatial integration in eq. ( |15f) requires some explanation since 
it does not appear in CKN: due to the finite thickness of the domain wall, the spatial region 
in which baryogenesis takes place is reduced. It was estimated in CKN that the baryon 
number violating interactions effectively turn off, in going from the symmetric to the broken 
phase, when the Higgs field reaches a value of g/An of its full VEV. For a profile of the form 
C (1 + tanh(z/A)) this occurs at z = zq = 1.5 A. If the Higgs potential is parametrized 
as \4> 2 {4> — 4> c ) 2 at the critical temperature, then A = (2/X) 1 ^ 2 (j)~ 1 , and the experimental 
lower limit on the mass of the Higgs boson gives V\ > 0.085. We assume that this bound 
is saturated, as well as the bound <p c > T c , which assures that sphaleron interactions in 
the broken phase are sufficiently suppressed to avoid washing out the baryon asymmetry on 
cosmlogical timescales. It is convenient to express zq in units of the characteristic width of 
the diffusion layer. The result is z^k^ = 1.0. 

In order to evaluate the density 5n of reflected quarks in front of the wall, one must 
assume a form for the reflection asymmetry AR(p z ). By analyzing existing results H0 it 
can be shown that AR is rather well fitted by the form 

AR(p z ) = N6(e)e a e- €/w ; e = p z - m. (16) 

The parameters N, w and a depend only on the quark mass, and a is neligible for a particle 
as heavy as the top quark [RJ. Then, although the explicit expression for 5n is rather 
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complicated, it takes a simple expression in the limit that the temperature is much larger 
than any other mass scales in the problem, 



s(zq) due to finite wall thickness is shown in figure 1. For the value of z determined above we 
see that s(z ) is approximately 0.12, which is a source of suppression that has not previously 
been taken into account. 

To compare to the Monte Carlo computation of 5n by CKN, we can use their deter- 
mination of the reflection asymmetry in our result, eq. fll7|) . For example, in the most 
favorable case when the width of the wall is assumed to be m _1 and CP violation is maxi- 
mal, they find a reflection asymmetry with N = 0.68 and w = 0.44m, from which we obtain 
Pb/s = 1.9 x 10 _7 /t(m/T c ) 4 . Our result is quite close to that of CKN for small wall velocities 
(v = 0.1) and m/T c < 1, but for larger m it begins to exceed their result. Of course in this 
region our approximations are no longer valid, but for a top quark of mass 170 GeV the 
Yukawa coupling is y = 1, and m/T c = ycf) c /T c = 1. 

We have computed the shape of the profile Sn(z) in the same approximation as used 
to obtain the integrated value of 5n above. It is shown in figure 2, where one sees that 
at a distance of 1 = 25/T it has dropped to 5% of its maximum value, and at 2/cq 1 to 
0.8%. Near the wall it falls off faster than an exponential, but asymptotically approaches 
exponential decay far from the wall. It is interesting to compare the width of the profile 
with the mean free path of quarks due to QCD scatterings, which we have computed to be 
A = 3T _1 . Figure 1 indicates that the profile falls by a factor of e in a similar distance, 
4T _1 , comparable to the diffusion length obtained using the position-space formulation of 
the diffusion problem |§. 

The behavior of 5n(z) we find is quite different from what would have been predicted 
from the diffusion equation, 




(17) 



after numerically evaluating the sum/integral appearing in eq. (0). The suppression factor 




(18) 
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where D z is the usual position-space diffusion coefficient of order 1/T, not to be confused 
with our D. The solution has the form 



n(z — vt) = n exp(— v(z — vt)/D z ) (19) 

for z > vt, the region in front of the wall. Thus in the small v limit the distribution 
approaches a constant in space, whereas ours has a finite extent (albeit vanishing height) 
even as v — > 0. This leads to at least one extra power of v~ l for the baryon asymmetry in 
the diffusion equation approach, due to the integration over z of the density in ©fl 

How do we understand the discrepancy between the two approaches? To bolster our 
confidence we can look for the same behavior in a simplified version of the same physical 
system, by working in one dimension and taking the infinite temperature limit, holding D 
fixed. Then eq. (|8]) becomes 

( % t #rl at = °> (2°) 



whose solutions are mere exponential and trignometric functions. In this case the integral 
for A~l analogous to (|13|) can be done explicitly and 8n(z) can be reduced to a single integral 
over u = (k/D) 1 / 2 , 

2N d r°° du _ zu 2 ( w + u 1 \ 

5n{z) = — / — e — , 21 

71 dw Jo u \w z + u A w + uj 

whose large distance behavior goes like z~ l l 2 . Again we see that the profile is not constant 

in position space even at zero velocity (although the length scale T 2 /D we had previously 

is now gone because of taking the T — > oo limit). However, we find that if we replaced the 

boundary condition flllD at the wall by one that did not depend on momentum, by letting 

the source term S(p z ) be a constant (let w — > 0), then the solution for would be a 

delta function at k = 0, and eq. fl2~T|) would indeed give a constant in z-space! So we have 

discovered one reason the two methods disagree: the Boltzmann equation demands that the 

boundary condition for the distributions at the bubble wall are functions of momentum, 

2 In fact there should be a second factor of w _1 because by Fick's law, the flux is J = —D z dn/dz = vuq 
at the bubble wall, so no is the flux at the wall divided by v. This extra factor does not appear in refs. Q 
or M. 
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whereas the diffusion equation has integrated over the momenta from the outset, and so 
cannot account for such details of the distribution. Of course the full Boltzmann equation 
should be the more accurate of the two methods. It would be interesting to compare the 
two predictions with the Monte Carlo method of CKN, but since the latter display data only 
for a single small value of v, it is impossible to measure their f-dependence in the region of 
v = 0.1. 

In summary, we have given an analytic expression for the baryon asymmetry due to 
quarks reflecting from domain walls during the first order electroweak phase transition, by 
solving the Boltzmann equation for the reflected particles. The solution is valid for slowly 
moving walls and temperatures which are large compared to the quark masses or the width in 
momentum space for the reflection asymmetry between the quarks and antiquarks. Although 
the top quark mass is actually not small compared to the critical temperature, the fact that 
we obtained quantitative agreement with the previous results of CKN for m/T as large 
as unity is encouraging since we don't expect it to be greater than this. Furthermore m 
vanishes (up to thermal loop contributions) in the symmetric phase outside the bubble, the 
region most important to the estimate of baryon production, so it might be hoped that the 
small-mass approximation is better than expected parametrically. An interesting feature of 
our expression is its insensitivity to the velocity of the expanding bubble wall, relative to 
predictions based on the diffusion equation. 

Our result makes it possible to estimate the baryon asymmetry AB directly in terms 
of the quark reflection asymmetry AlZ(p z ), an independent calculation which is being re- 
peated || because of apparently conflicting calculations in the literature @0. The quan- 
titative results of electroweak baryogenesis through quark reflection hinge upon this since 
the momentum-space width of AlZ(p z ), on which AB depends linearly, falls quite sharply 
with increasing quark mass, whereas for small masses AB is suppressed by a factor of m 3 . 
Further application of the present work will be given in ||. 

I would like to thank Larry McLerran and Sonia Paban for important contributions 
during the early stage of this work, which was supported in part by DOE grant DE-AC02- 
83ER-40105. I also thank Kimmo Kainulainen and Alejandro Ayala for helpful discussions. 
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Fig. 1: logarithm of the suppression 
factor due to finite wall thickness, zq is 
proportional to the bubble wall thickness 
(see text), and is in units of the diffusion 
length A^ 1 = 25T- 1 . 
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